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Abstract: In this paper, we put constraints on neutrino properties such as mass and 
degeneracy parameters from WMAP5 data and Hght element abundances by using a 
Markov chain Monte Carlo (MCMC) approach. In order to take consistently into account 
the effects of the degeneracy parameters, we run the Big Bang Nucleosynthesis code for 
each value of and the other cosmological parameters to estimate the Helium abundance, 
which is then used to calculate CMB anisotropy spectra instead of treating it as a free 
parameter. We find that the constraint on rui, is fairly robust and does not vary very much 
even if the lepton asymmetry is allowed, and is given by ^ tjIu < 1.3 eV (95%C.L.). 
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1. Introduction 

Neutrino masses are tlie key feature beyond the standard model of particle physics because 
neutrinos are assumed to be massless in the minimal standard model. However, oscillation 
experiments recently suggest tiny but non-zero masses of neutrinos. Unfortunately they 
can probe the mass (squared) differences between different mass eigenstates but cannot 
determine absolute neutrino masses. In addition, the mass of electron-type neutrino is 
bounded as rrijy^ < 2.0 eV (95% C.L.) from the current laboratory experiment of tritium 
beta-decay Q. However, the most stringent constraint on the sum of neutrino masses comes 
from cosmological observations. For example, by using the SDSS luminous red galaxies Q, 
the sum of neutrino masses is constrained as '^nii, < 0.9 eV (95% C.L.), though there 
is still uncertainty in the galaxy bias. On the other hand, Ichikawa et al.|3|, ^ showed 
that similar constraint can be obtained by CMB data alone, which is free from such an 
uncertainty and the most conservative. The latest WMAP results yield X^mjy < 1.3 eV 
(95% C.L.) for the ACDM model [|[. This constraint is shown to be robust over the 
different cosmological models, even if we abandon the assumption of the spatial flatness 
of the universe, change the equation of state of dark energy, and include tensor modes 

Another assumption in the standard cosmology that is usually taken is the symmetry 
in the lepton sector. This is closely related to cosmological neutrinos since neutrinos have 



all the lepton number of the universe due to its electronic charge neutrality. The standard 
sphaleron process forces the amount of lepton asymmetry to be the same order as that 
of baryon asymmetry. However, several ways are proposed to avoid such sphaleron effects 
1^, ^, ^, IC, 11, 12, 1^, 14, 15, !(:]. For example, if Q-balls are formed after the Affleck-Dine 



leptogenesis, the large lepton asymmetry inside Q-balls is protected from sphaleron effects 
and the unprotected small lepton asymmetry is converted into the small baryon asymmetry 



12]. Thus, small baryon asymmetry is naturally compatible with large lepton asymmetry. 



In addition, observational constraint on lepton asymmetry is much weaker than that of 
baryon asymmetry. Then, one may wonder how robust the constraint on neutrino masses 
is in the possible presence of large lepton asymmetry. This is the main topic we would like 
to pursue in this paper. 

The lepton asymmetry which is usually denoted by the degenerate parameter ^ is an 



interesting quantity to be constrained |17, 18, 19, [2^, where ^ is defined as a chemical 
potential between neutrino and anti-neutrino normalized by the neutrino temperature. 
Recently, Popa and Vasile have obtained a limit on ^ by using the WMAP 5 year data 
combined with the other CMB data and the LSS data ||2^. However, since there are 
uncertainties of biasing and non-linearity as regards galaxy clustering data, it would be 
useful to constrain ^ from WMAP5 alone. We also consider indirect effects of ^ on CMB 
through the primordial helium abundance Yp generated during the process of the Big Bang 
Nucleosynthesis (BBN). 

In this analysis, we consider three light neutrinos, based on the result of high energy 
experiments such as LEP, i.e. N,y = 2.984 it 0.008 |^l|]. We assume equal neutrino masses 
between fiavors, such as = m^^ = rriy^ = rriy^, which yields ^ = 3m^. It is because 
the accuracy of the present cosmological observation is not yet up to specify their mass 
differences. Potentially, the degeneracy parameter ^ can take different values for each flavor 
of neutrinos. In this paper we consider the two cases: one is the case that = Ca* = ^r, 
and the other is that = = / ^g- The flrst case is motivated by Refs. |22, 22, 24 1, 
which show that if a flavor mixing angle of neutrino ^13 is large enough, each degeneracy 
parameter becomes same before BBN, namely = (,fi,r- However, as recently stressed 
by Ref. |25|, while this conclusion is basically correct, oscillations take place so late that 
perfect equilibrium is never achieved and the ratio of imperfectness strongly depends on 
the magnitude of ^13 [25|. It is also pointed out that if a Nambu-Goldstone particle such as 
majoron on interacts with neutrinos, flavor mixing can be suppressed and leads 7^ Cm.t 
1 26]. Therefore, we also consider the second case 7^ C^.r- The second case is also 
motivated by the fact that and have the same effects both on the BBN and CMB 
because we consider the degenerate masses, while does not. 

This paper is organized as follows. In the next section, we summarize effects of mu and 
on BBN and CMB. In Sec. ^, we give the cosmological model that we use in parameter 
searching of MCMC analysis. In Sec. we present our limits from the observational data 
based on WMAP5 and explain the correlation between the bound of each parameter. In 
final section ^, we provide the summary of this paper. Throughout this paper, we set 
c = h=l. 



2. cosmological effects of and m, 



In this section, we review the effects of the neutrino mass rrii, and the degeneracy parameter 
on the predicted hehum abundance Yp in BBN theory and the predicted CMB power 
spectrum Ci. The degeneracy parameter of ith. flavor of neutrino is defined as 

C^ ^ ^, (2.1) 

where is the chemical potential. 
2.1 Effects on BBN 

The observed primordial light-element abundances constrain the conditions during the BBN 
epoch from the time of weak reaction freeze out {t ~ Isec, T ~ IMeV) to the freeze out of 
nuclear reactions {t ~ lO'^sec, T ~ lOkeV). In the standard BBN theory, the abundance 
of light elements is determined by only one parameter, namely, baryon-to-photon ratio. 
However, the lepton asymmetry changes the abundance in two ways. 

One is that the extra energy density coming from the neutrino degeneracy changes the 
expansion rate of the universe during the BBN epoch and hence the primordial abundances. 
Specifically, the Fermi-Dirac distribution functions of neutrinos and anti-neutrinos with 
neutrino degeneracy are given by 

A(g,^») = i^e^-^. ^ M{Q,Ci) = (2-2) 

where q is comoving momentum written by proper momentum P and neutrino temperature 
Tjy as g = P/Tiy, and is the degeneracy parameter of ith flavor of neutrino. The energy 
density for massless degenerate neutrinos is given by 

Pu + Pp = ^(/cfiT^)^ / + fuM^ii)) = ^{Puo + Puo) > (2-3) 

i 

where p^/g is the energy density for massless non-degenerate neutrinos and N^q is the 
effective number of relativistic degrees of freedom of neutrinos given by 



In the non-degenerate neutrino case, N^fi = 3. Therefore the neutrino degeneracy increases 
the energy density of neutrino. 

The other is that positive suppresses p — >■ n weak reactions and leads the neutron- 
to-proton ratio to smaller value both in and out of the chemical equilibrium. Therefore 
the abundance of helium Yp decreases monotonically with increasing ,^e- Note that the 
dependence on mjy can be neglected because rriu is much smaller than the temperature at 
the BBN epoch and does not affect the reaction rate. 



We used the Kawano BBN code [27| to include the two effects described above and 
to calculate the abundance of light elements correctly. We also updated the nuclear re- 
action rates [28|. In Fig. |l[ the effects of and ^^^^ on the helium abundance Yp are 




Figure 1: Effects of the degeneracy parameters, and Cm-t' '^^ the hehum abundance Yp. Here 
we fix the baryon abundance uji, to the mean value from WMAP5 alone analysis for a power-law 
fiat ACDM model [|. 

demonstrated. The larger A'cfr makes the universe expand faster and more neutrons sur- 
vive before the nucleosynthesis begins. Therefore the helium abundance becomes larger as 
^^^r increases. On the other hand, Yp decreases monotonically as increases because the 
chemical potential of electron type neutrinos suppresses the neutron-proton ratio at the 
onset of BBN as mentioned above. 

2.2 Effects on CMB 

Cosmological Microwave Background (CMB) photons are the remnants from the Big-Bang 
in the early universe. They have evolved through the interaction with the matters and 
the gravitational potentials, and have been released from the last scattering surface at 
cosmological recombination (z ~ 1088) to reach us through the large scale structure of the 
universe. Therefore CMB contains plenty of information about the energy contents of the 
universe. In the simplest cosmological model, the universe is assumed to be spatially flat, 
and neutrinos are massless and not degenerate. The angular power spectrum of the CMB, 
Ci, can be described by the six cosmological parameters, namely, the energy density of 
baryon and dark matter, the present Hubble parameter, the amplitude of the primordial 
curvature perturbation and its spectral index, and the optical depth due to re-ionization. 
However, if rrijy and are finite values, both of them also affect Ci both at the background 
and the perturbation levels. 

In the background evolution, the mass of neutrino tjIu and the neutrino degeneracy 
determine the energy density and pressure of neutrinos. The energy density of massive 
degenerate neutrinos is described as Pu + pp = '^i'l^vi'^'Uii^i) + i^uiid)), where n^. is the 
number density of neutrinos. Therefore the energy density parameter of them, which is 
the normalized energy density by the critical density pcrit = (3i^^)/(87rG), is given as 



o 
o 

O 

o 
<x> 

o 
o 
o 



/Ocrit(O) 



(2.5) 



-4- 



Here the argument (0) means the value at present and the nth moment is defined as 



&n{C^) q^dqiU^ + M)/ . (2.6) 

When we can regard neutrinos as massless degenerate particles, pu + pp is expressed by 
using just like the discussion in the previous subsection. 

Meanwhile, in the perturbation evolution, the changes on the gravitational source 
term of the coUisionless Boltzmann equation of neutrino. Specifically, its coefficient, df /dq, 
is modified as 

rfln(/i., + fpj) ^ q 1 + coshgjcoshg 

ding e~'> + cosh,^i cosh^j + coshg 

Because we consider the situation that neutrinos have degenerate masses, and because neu- 
trinos are coUisionless at the CMB epoch irrelevant to their flavors, there are no difference 
in the effects on Q between and S^^^r contrary to the product Yp at BBN. For the details 
of the effects of rriu and ^ on CMB we refer readers to the references p9| , |30t| . Here we 
derived above equations following their notations. This dependence of the perturbation 
on ^ is absent for the massless degenerate case, because the momentum dependence in 
distribution function is integrated away into the energy density. In this case, the effect of 
on Ci is identical to that of N^ff through Eq. (p.4|). Although this momentum depen- 
dence should have some distinctive effects on the evolution of perturbations in the massive 
degenerate case, we see that the dominant effect still comes from the change in A'cfr for the 
neutrino masses X^rrij^ < 1.5 eV considered in this paper. 

We inserted these two modifications in the Boltzmann Code for Anisotropies in the 



Microwave Background (CAMB) pl|, 32, 33]. Then we checked the consistency between 



the Ci's given by and those calculated using our modified code. In Fig. effects of 
nil, and ^ on the CMB temperature power spectrum Cf^ are demonstrated. 

Now we explain the dependence of the CMB spectrum on and m^. If the effective 
number of neutrino, N^s, increases by increasing the sound horizon at the last scatter- 
ing becomes smaller and the epoch of radiation-matter equality delays leading the larger 
decay of the gravitational potential. These effects shift Ci to the higher / and derive the 
enhancement of the 1st peak. In addition, through the effects of free streaming of these 
ultra-relativistic neutrinos, the smaller scale gravitational potentials have been damped 
just prior to the beginning of the acoustic oscillations of the baryon-photon fiuid. This 
causes the smaller temperature fluctuations. These effects on Ci from N^s are discussed in 
detafl in pi. 



Next we focus our attention on the dependence on rrii, In massive non- 

degenerate case (^j = 0), neutrino mass of 0(eV) causes overall horizontal shift and sup- 
pression around the first peak (when we fix ivh + and change uj^ to keep flatness). The 
horizontal shift comes from the fact that the larger niu (more non-relativistic particles at 
present epoch) implies that the distance to the last scattering surface is shorter and the 
peaks move to smaller /. However, this shift is mostly canceled by the downward shift 
in Hq. Therefore this does not produce a neutrino mass signal. If rriu > 0.6 eV, massive 
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Figure 3: Effects of ^ on Cf^ for the case = 



Effects of and ^ on Cf^ for tlie ^^^r witli massless neutrinos. Two dotted lines 
,r with Yp fixed as 0.24. Here we fix and a dashed line are Cf^^s affected by Yp related 

with ^ at BBN. Solid lines show Cf^'s affected by 
Yp unrelated with ^. Here we fix uJb + tOc- Other 
cosmological parameters are taken as the mean 
value from WMAP5 alone analysis for a power- 
law fiat ACDM model |] i.e. + f^^ + f^A = 1- 



Figure 2: 

case = 

ujb+ujc and change oju- Other cosmological param- 
eters are taken as the mean value from WMAP5 
alone analysis for a power-law flat ACDM model 
II i.e. f7t + f7c + + f^A = 1. 



neutrinos on average become non-relativistic before the epoch of recombination and only 
in this case, the neutrino mass can significantly imprint a characteristic signal in acoustic 
peaks (specifically, the matter-radiation equality occurs earlier due to less relativistic de- 
grees of freedom and 1st peak is suppressed by the smaller early-integrated Sachs- Wolfe 
effect) |, ||. 

In massive degenerate case (^j ^ 0), the behavior of C; can be understood by combining 
the effects in the massless degenerate case and the massive non-degenerate case. In Fig. ^ 
we depict the CMB angular power spectrum to illustrate the effects of mass and degeneracy 
of neutrinos. We see that the m,y-dependence does not much depend on ^. As mentioned 
above, roughly speaking, neutrino mass constraint from (WMAP-level) CMB observation 
is determined by a critical mass 'irriy^c which is defined by z^ec — Znr, where Znr denotes 
the redshift when neutrinos on average become non-relativistic. When there is a finite 
the neutrino momentum distribution is modified and so is this critical mass. In this case, 
the critical mass is obtained by 

"^-^^^ - ^^^^^ ^F^O, (2.8) 



m,,,(e = 0) P,(C = 0) 

where Pu{£,) and PuiC = 0) are the average momenta of neutrinos in massive degenerate 
and massive non-degenerate cases, given by P^, ^ Pu/^u, and Fi, is given by 





Figure 4: The ratio of the neutrino masses between degenerate and non-degenerate cases for 
which the neutrinos become non-relativistic at recombination. 



Numerically we depict Fi, in Fig. ^ From the figure we found F^, < 1.14 for < 2. 
Therefore the critical mass {3m,y^ciS, = 0) = 1.7 eV for massive non-degenerate case) shifts 
to, at most, 1.9 eV for massive degenerate case. As it is indicated Q that niiy is already 
limited less than m,y^c, we can expect that the upper bound of ^ from CMB is almost 
invariant even in the lepton asymmetric universe, as will be shown later. 

Finally, we briefly explain the dependence of the CMB power spectrum on Yp. The 
main effect of Yp on the CMB power spectrum comes from the diffusion damping at small 
scales. When Yp is larger, the number of free electron becomes smaller, because the max- 
imum number of free electron in the fully ionized universe is propotional to (1 — Yp/2). 
The number of free electron becomes even smaller at the epoch of hydrogen recombination, 
since it is easier for electrons to recombine with ^He than with H and consequently more 
electrons have been captured into the ^He nuclei by that epoch. Thus the Compton mean 
free path becomes larger for larger Yp, which mean that the diffusion length of photon 
becomes also larger. Since the photon-baryon tight coupling breaks down at the photon 
diffusion scales, the fluctuation of photon is exponentially damped due to the diffusive mix- 
ing and re-scattering. Hence the CMB power spectrum is damped more for larger values 



oiYp. Detailed explanations are found in Refs. |36, 37, 38]. This effect is rather small and 
not relevant to WMAP5 level observation. However, in our case, since Yp can take very 
large or small value with ^ ~ 0(1), this may affect constraints. For example, two models 
with = C/i,r = 1-5 and —1.5 give the same CMB spectra if we fix = 0.24. However, 
if we consider the effects of ^ at BBN epoch correctly, the two models predict Yp = 0.048 
and 0.780, respectively as shown in Fig. |l]. Therefore, the model with negative degeneracy 
parameter should lead the larger diffusion damping as shown in Fig. ^. The CMB power 
spectrum can be affected by ^ through such an indirect manner [39|. 



3. Model and likelihood analysis 



In what follows we will put constraints on the neutrino degeneracy parameters and 



masses along with the other standard cosmological parameters. For this purpose we 
implement CAMB and BBN codes modified by us for lepton asymmetric cosmology in the 
CosmoMC Monte Carlo Markov Chain (MCMC) public package [^]. In our analysis we 
use observational data set in several combinations; WMAP5 alone, WMAP5 + prior on 
Hubble parameter (Hq = 72 it 8 km/sec/Mpc (68%C.L.)) given by Hubble Space Telescope 
(HST) and WMAP5 + "^He and D/H abundances (BBN). Throughout our analysis 
we assume that the universe is spatially flat and neutrinos are degenerate in mass. We 
explore the likelihood in 10 dimensional parameter space, i.e., 

p= {LUb,UJdm,A,ns,9s,T,Asz,^e,^ti,r,'^u)- (3.1) 

Here u!/,, uj^m and lUi, are the energy densities of baryon, dark matter, and massive neutrino, 
respectively, Ag and Ug are the amplitude and the spectral index of the primordial density 
perturbation power spectrum at the pivot scale /c* = 0.002 Mpc~^, Og is the acoustic 
peak scale, r is the optical depth of reionization, and A^z is the amplitude parameter 
of thermal Sunyaev-Zel'dovich (SZ) effect |^^. Note that uj^, is included in ujdm so that 
uJdm = oJc + ^v, where ujc is the energy density of cold dark matter. To calculate Ci in the 
parameter estimation, we estimated Hq from Og by using the fitting fomula in ACDM for the 
redshift of recombination Zrec- This is the default method in CosmoMC. In the degenerate 
case, we should consider the effect of ^ on this fitting formula. However, as we found that 
the shift of Zrec for ^ = 2 is less than 0.1% from the result of the fitting fomula in ACDM 
for ^ = 0, we used it safely even in the degenerate case. As for the helium abundance Yp 
used in CMB calculations, we consider two situations: one is that we (approximately) fix 
Yp = 0.24 regardless of the cosmological parameters, and the other is that Yp is derived at 
each time of MCMC step from the other cosmological parameters using the modified BBN 
code. Hereafter we denote for the latter case the helium abundance as 1^^^^, i.e., 

yp = y/BN(a;fe,ee,W) • (3-2) 

Note that for the former case we vary the degeneracy parameter ^ freely in the MCMC 
analysis even if the helium abundance is kept fixed as Yp = 0.24. 

In the usual CMB constraints the helium abundance is often fixed to the standard value, 
say Yp = 0.24. This condition is justified by the fact that current CMB measurements have 
already put a tight constraint on w^, and for this constrained range of the CMB spectrum 
does not change very much even if we allow the helium abandance to vary according to the 
standard BBN prediction. 

However, in the lepton asymmetric universe considered in this paper, the value of Yp 
can deviate significantly from Yp = 0.24 depending on ^ as discussed above. In this sense, 
any analysis of CMB constraint on ^ with fixed Yp is inconsistent unless ^ is constrained 
tightly. In our CMB power spectrum calculation as is used in later analysis, Yp is set by the 
BBN theoretical calculation as a function of w;,, and C/i,T- We also investigate constraints 
in the lepton asymmetric models with fixed Yp for comparison. 



4. Results and comparisons 



In this section, we describe our results and give their interpretations. For a test calculation 
we first ran CosmoMC to constrain neutrino masses from WMAP5 alone with = 0. We 
obtained consistent but slightly tighter result, J2^v < 1-2 eV (95%CL), while WMAP 
team obtained 1.3 eV p. This difference would be attributed to differences in the treat- 
ment of cosmological parameters as well as their priors |44]. In general, posterior 



distributions depend on prior distributions in Bayesian estimation pO| ]. 

Next, we ran a MCMC with all of the parameters to be fixed other than ^^^^ and 
confirmed that the ID posterior distributions of ^^^^ are symmetric about ^ = 0. Using 
this fact we restricted the prior range of ^^^r only to positive values, i.e., we actually 
estimate \(,^^t\ from MCMC analysis. 

We perform MCMC analyses for three combinations of observational data sets de- 
scribed below. In all analyses, we neglect effects of CMB lensing because they slightly 
change Ci but do not change the posterior distributions very much. 

4.1 Constraints from WMAP5 alone 

We first present the result of constraints on cosmological parameters from WMAP5 alone 
in TABLE ||. In that table, we focus on 5 parameters, ^j, m,y, Hq, Yp, uJdm- For cosmological 
models we consider the following four cases and give the results separately. We present the 
mean values and 68% and 95% confidence regions for the case .^e = ^^^,T in the left panel 
of TABLE 0, and for the case 7^ ^fi,T in the right panel. For each panel we separate the 
results for cases whether we use BBN code to calculate Yp or we approximately fixed Yp 
value to 0.24. The results for the massless neutrino cases are given in the appendix. 

As is found in TABLE the constraints on ^ from CMB are highly dependent on 
whether one fixes helium abundance or derives it from BBN theory. First, let us consider 
the case with fixed helium abundance. In this case, the constraints on ^ from CMB becomes 
symmetric about = 0, because positive and negative ^ give identical effects on the CMB 
spectrum through Eqs. (|2.4| ) ( |2.5D , and ( |2.7D . In reality, however, the positive and negative 
^ give different effects on the CMB spectrum through the helium abundance. In CMB 
theory, helium abundance is related to the number density of free electrons, and hence the 
Silk damping scale depends on it. So the CMB can put constraint on Yp, and therefore 
is further constrained through this effect if one derives helium abundance from BBN 
theory for the CMB calculations. From the fact that the current upper limit on helium 
abundance from CMB is found to be < 0.44 and there is no lower bound |38|, negative 



is disfavored, because negative leads to large Yp as is discussed in As for the 
other cosmological parameters, such as ujh and n^, are tightly constrained for all cases 
and the constraints do not change between different models considered here, because the 
dependence of Ci on these parameters does not degenerate with that on and Yp. 

One may find it inconsistent that the constraints on and ^^^^ with fixed Yp = 0.24 are 
different because in this case and ^^^^ play the same role in the CMB. However, in fact, the 
bound of ^^^r should be tighter than that of because we implicitly assume the condition 
^^l = in this analysis. We analytically estimated this difference as ICc/Cai.tI ~ 1-3, 











parameter 


Yp = 0.24 


V — vBBN 
— -fp 


Yp = 0.24 = Y^^^ 



lOOwfe 2.205 ±0.070 2.206l|]-[574 2.2mtQ^l 2.199 ±0.075 

n 1 C9Q+0.015 1 7QQ+0'.0079 n 91 1 1 +0-0202 p, r,-, cro+0.0208 

< 2 41 1 20+°-5°+^-39 < S 47 1 fi^+0.63+2.00 

: -_ < 2.70 < 2.81 

Xl"!^ (eV) < 1.2 < 1.3 < 1.2 < 1.3 

Yp - < 0.36 - < 0.32 

Table 1: Mean values and 68% C.L. errors of the cosmological parameters obtained from the 
analysis of WMAP5 alone in the lepton asymmetric universe with massive neutrinos. For the C — | 
constraints on ^, 95% C.L. errors are also shown, and the upper bounds on ^, ^ m^, and Yp are at 
95% C.L. The left side of the table is in the case of = Cm,^- "^^e right side is in the case 7^ £,^^,T 
Two parameters at the bottom are the derived parameters from the MCMC sampling. 



o 

which can be derived as follows. If we varied each of the degeneracy parameter freely, 
the distributions for in and should have been identical in the case where helium 
abundance is kept fixed and the masses are degenerated. The effects of ^ in this case 



mainly come through an increase of neutrino energy density, i.e., Eq. (2.4). Therefore, 
once the constraints on can be obtained in the form as 



o 

with some value of a, then, the constraints on should read dZ) 

o 



30 /^WV , 1^ /^^'' 
7Vvr/ 7Vvr 



< a . 



(4.2) O 



We found that this estimate is in good agreement with our MCMC constraints. The result 
of the constraints from WMAP5 alone is depicted on Hq-Y^ rrii, plane in Fig. ^. In that 
figure, the points are drawn such that the density of the points is proportional to the 
number of the MCMC samples (i.e., probability distributions), and colors represent the 
value of ^. One can see from this figure a well-known Hq-^ degeneracy with fixed 
^ values. We find that the upper bound of is not changed significantly because 

this upper bound is lower than the critical mass 2>mu,c- It is interesting to note that the 
massive neutrino cosmology with ^ rriy ~ leV is compatible with the Hubble parameter 
Hq 72 km/sec/Mpc if one allows ^ 1 lepton asymmetry, as will be discussed in the 
next subsection. 

From WMAP5 alone, the upper bounds of are actually determined by the flat prior 
on the Hq which is taken as Hq < 100 km/sec/Mpc. Therefore, analysis which takes into 
account external prior on Hq should be performed. In the next subsection, we apply a 
prior of Hq from HST. 
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Figure 5: The constraints in the 3D region of ^TOi^, Hq and ^ from WMAP5 alone in the case 
of = with Yp = 0.24 (the left figure) or Yp = Y^^^ (the right one). The density of the 
points, whose colors are coded according to the value of ^, represent the probability distribution 
from MCMC samples. The contours are the 68% and 95% confidence regions. 



4.2 Constraints from WMAP5 + HST 

In this subsection, we present the constraints from WMAP5 + HST, with the helium 
abundance beinff set to Y^^^ (see Eq. (jS^)). Here the prior from HST is Hq = 72 ± 8 
km/sec/Mpc (68%C.L.) @. These results are presented in TABLE |. 

In TABLE |^, we found that all the parameters except for ^ rriiy and Yp are limited 
more tightly than those limited from WMAP5 alone. This is due to breaking the Hq - ^ 
degeneracy. 

In the standard non-degenerate case, by adding HST prior the upper bound of ^ rrii, 
becomes tighter than that from WMAP5 alone, and we found ^rriu < 1.0 eV. This 
is because the prior of HST breaks effectively the well-known ^ - Hq degeneracy. 
Interestingly, however, we found that the upper bounds of ^ m,^ from WMAP5 -|- HST do 
not change very much compared with the upper bounds from WMAP5 alone if the lepton 
asymmetry is allowed. In the degenerate case, higher neutrino masses are compatible with 
the Hubble parameter around Hq ~ 72 km/sec/Mpc with non-zero as depicted in Fig. 
??. This is because the shift of the peak to the lower / due to larger rrii, can be compensated 
by the shift due to larger ^, with the Hubble parameter fixed. Therefore, HST prior does 
not improve the limits on neutrino masses in the lepton asymmetric universe. 



4.3 Constraints from WMAP5 -|- BBN 

Finally, we present the constraints from WMAP5 combined with observations of light 
element abundances of 1^ and D/H. We use only Yp = 0.250±0.004 |Q in the case = Caj.t- 
In the case 7^ ^^l,T^ there is a parameter degeneracy between and ^^^^ for the constraint 
from the helium abundance alone. Therefore, we further use the deuterium abundance 
constraint as D/H = (2.82 ± 0.27) x 10~^ [46] in this case to break the degeneracy. The 
results are listed in TABLE. ^, and our constraints on the sum of neutrino masses are 
1.2 eV and 1.1 eV respectively for the cases = Caj.t and 7^ C/i,r- They do not depend 



parameter C = Ce = ifi,T Ce + 



lOOcJb 2.202t°-°^^ 2.193t[![]^^ 

niR91 +0-0045 -I q 4 7+0'.0132 

c n Qfi+o.50+1.35 1 40+0.55+1.92 

?e '-'■^"-0.40-1.52 -'-•^°-0.64-1.68 

IWI - <2.50 



(eV) < 1.3 < 1.4 

< 0.41 < 0.32 

Hq (km/sec/Mpc) 71.8+3;^ 74.0+3^ 

Table 2: Mean values and 68% C.L. errors of the cosmological parameters obtained from the 
analysis of WMAP5 + HST in the lepton asymmetric universe with massive neutrinos. For the 
constraints on ^, 95% C.L. errors are also shown, and the upper bounds on ^, ^m^, and Yp are 
at 95% C.L. The left side of the table is in the case of — f^.r. The right side is in the case of 
7^ C^i.T- Here we use the relation Yp = Y^^^ in the calculation of the CMB power spectra. Three 
parameters at the bottom are the derived parameters from the MCMC sampling. 




Figure 6: The constraint in the 3D region of ^ : 
one) from WMAP5 + HST in the massive ^ 5,^, 
are depicted in the same way as Fig. 5. 




<.uj Ho and (the left figure) or ^^^t- (the right 
case with Yp = Y?^^. The points and contours 



on the assumption of = ^^l,T and 7^ ^m.t- Again, we have found that neutrino mass 
constraints do not change much depending on the assumptions on ^. 

As for the constraints on the degeneracy parameters, we obtain ^ = — 0.013+g Qj^7 in the 
case ie = C^l,T, and = 0.034lo;o25 and \C^i,t\ < 1-60 in the case ^ C/^.r from WMAP5 + 
BBN. In relation to them, we also performed a conventional analysis. Namely, we searched 
three-dimensional parameter space of ojb, and ^^^t- by BBN calculation. Here, we take 
the baryon density of uJb = 0.0225 ± 0.0006 obtained from WMAP5 data in flat ACDM 
model. The constraints, ^ = -O.OlSto;!]^^ in the case = C^.r, and = 0.002t[J;[J^^ and 
|^/i,r| < 1-61 in the case 7^ C/^.r are consistent with the cases above confirming the validity 
of the conventional procedure. 
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parameter ^ = = g/.,r / ^^l,T 



lOOcJb 2.215 ±0.066 2.221 ±0.061 

n 1 1/1 7+0- 0079 n 1 -^i 1+0-0056 

^dm -0.0078 '-'■-'-•-'-'^-'^-0.0096 

(: _n ni'^+0 016 -t n 033 O nQ4+0.017+0.075 

IWI - < 1.60 



(eV) < 1.2 < 1.1 

0.250 ± 0.004 0.250 ± 0.004 
Hp (km/sec/Mpc) 66.4+^;^ 70.0^3;^ 

Table 3: Mean values and 68% C.L. errors of the cosmological parameters obtained from the 
analysis of WMAP5 + BBN in the lepton asymmetric universe with massive neutrinos. For the 
constraints on ^, 95% C.L. errors are also shown, and the upper bounds on ^ and 'Y^rrii, are at 
95% C.L.. The left side of the table is in the case of = ^tJ.,T- The right side is in the case of 
7^ ^/i.T- Here we use the relation Yp = 1^^^^ in the calculation of the CMB power spectra. Three 
parameters at the bottom are the derived parameters from the MCMC sampling. 



5. Summary 

In this paper, we have discussed the effects of neutrino masses and neutrino degeneracies 
on CMB and BBN, and investigated how the latest observations put constraints on these 
properties simultaneously. In particular, we have examined how robust the constraint on 
neutrino masses is in the lepton asymmetric universe. Our constraints on the neutrino 
masses are rriy < 1.3 eV at 95% C.L. from WMAP alone. The constraint does not change 
very much even if we put the HST prior in the lepton asymmetric universe, as explained in 
Sec. 4.2. We investigate the two types of lepton asymmetry, i.e. the cases with = ^fi,T 
and 7^ '?/i,r5 and found that the constraints are almost the same for both cases. We also 
have taken into account the measurements of light element abundances, and found that 
changes in constraints are very small. These results are summarized in the left panel of 
Fig. |. 

As for constraints on the degeneracy parameters, we have obtained at 95% C.L. 
—0.41 < < 2.59 for = ^fM,T- In the case of .^e 7^ C^J.,T, the constraints are —0.19 < 
< 3.65 and [^^,t1 < 2.81 from WMAP5 alone. We point out that they are limits ob- 
tained by imposing the BBN relation for Yp and this helps to place tight lower bound on 
^e- A part of these results are shown in the right panel of Fig. 0. There, we can notice 
that the posteriors are skewed. In addition, the right panel indicates that ^e's prefer the 
positive value, however, they are still consistent with = 0. Although these constraints 
are much weaker than those obtained from Yp measurements in HII regions and direct use 
of BBN theory, this would be more useful in future CMB measurements such as Planck 
||39| , poll to constrain ^g- 
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A. Constraints for the massless degenerate case 

In this Appendix, we present the result of constraints on cosmological parameters for 
massless degenerate case with 1^'s relation as 1^ = 1^^^^. In TABLE ^ and |5|, we show 
all the results. In this case, basically, the bounds are determined by the similar effects 
as we discussed in the massive case. Here we stress the point simplified by dealing with 
neutrinos as massless particles is that Hq is completely determined by the degeneracy effect 
of Hq — reported in ||3^ once ^'s is limited by the effect of Yp. 
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